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1—1 Abstract 

Qh Imaging techniques give a fundamental support to medical diagnostics during 

Qh the pathology discovery as well as for the characterization of biological struc- 

S tures. The imaging methods involve electromagnetic waves in a frequency 

^ range that spans from some Hz to GHz and over. Most of these methods 

• involve scanning of wide human body areas even if only small areas need 

O to be analyzed. In this paper, a numerical method to evaluate the shape 

^ of micro-structures for application in the medical field, with a very low in- 

r-| vasiveness for the human body, is proposed. A flexible thin-wire antenna 

^ Ph radiates the VHF waves and then, by measuring the spatial magnetic field 

distribution it is possible to reconstruct the micro-structures image by es- 

^ timating the location of the antenna against a sensors panel. The typical 

,-H inverse problem described above is solved numerically, and first simulation 

^^ results are presented in order to show the validity and the robustness of the 

(^ proposed approach. 
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1. Introduction 

The support given by the imaging to medical diagnostics is fundamental 
during the pathology discovery as well as for biochemical characterization 
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of biological structures PP-[Z]. The imaging methods involve electromag- 
netic waves in a frequency range that spans from some Hz to GHz and over. 
Hence, the understanding of biological structures response to the electro- 
magnetic field is fundamental. The investigation of the interaction between 
healthy or pathological biological tissues and electromagnetic waves is a hot 
topic yet [S]. In fact, understanding how an electromagnetic wave interacts 
with the human body becomes increasingly important if one considers that 
the most of imaging methods involve scanning of wide areas of the human 
body, even if only small areas need to be analyzed. In this way, although a 
wide scanning allows to acquire a big amount of data by single exposition, 
also areas of the body not interested in the diagnostics are exposed to the 
waves, so increasing the invasiveness for the patient. For this reason, new 
imaging systems able to analyze only the area under diagnostics (confined 
scanning), are becoming more and more investigated. A big input to micro- 
imaging systems has been given by the micro-electronic technology, which 
allows the development of systems for the scanning of confined small areas. 
In fact, emitters and receivers (sensors) smaller and smaller can be made, so 
that the final size of the imaging systems as well as the waves spot become 
very small. Generally, the acquired data need to be appropriately elaborated 
extracting the imaging information by means of inverse optimization algo- 
rithms |3],|1],[Z!- Through these algorithms high quality information can be 
extracted from the data acquired by the sensors, even if the quality of the 
sensors signal is low, because of their small area which may determine a sig- 
nificant amount of noise. In this paper a method to elaborate the shape of 
microstructures for application in the medical field is proposed. The method 
works with low power waves in the Very High Frequency (VHF) range, in 
order to achieve a low invasiveness for the human body. The method uses a 
system endowed with: a microtransmitter to emit a magnetic field, a sensors 
panel to acquire the spatial distribution of the magnetic field and an elabora- 
tion logic to acquire and elaborate the sensors signals. The microtransmitter 
radiates the VHF waves by means of a microantenna, which is able to take 
the shape of the target structure. If the micro-antenna is assumed to be 
a sequence of thin-wire interconnected short dipoles, then it is possible to 
reconstruct the micro-structures image by measuring the spatial distribution 
of the magnetic field. In fact, the shape reconstruction is possible by es- 
timating the location of thin-wire antenna against the sensors panel. The 
recognition problem of the thin-wire antennas location by magnetic field can 
be addressed as an inverse problem. The thin-wire antenna is supposed to 



be a sequence of linear segments: given a model for the characterization of 
the magnetic field at a point in space (forward problem) and given a set of 
measurements of the magnetic field amplitude, it is possible to solve the in- 
verse problem in terms of the distance of the antenna from the sensors panel. 
In this paper, preliminary results about the previous basic idea, are reported 
concerning two simulated scenarios. Namely, the spatial distribution of the 
emitted magnetic field is simulated through a numerical model based on the 
method of moments (MoM), where the first kind Fredhloms integral equation 
is solved by the point matching procedure [9]-[T5]. The spatial distribution 
of the magnetic field is evaluated on an area equivalent to the area of the 
sensors panel. These magnetic field values are used as the measured input 
for the inverse problem. The Levenberg-Marquart algorithm [16] , [17] is con- 
sidered in solving the inverse problem by means of the minimization of the 
Euclidean distance between the measured field and the field generated by a 
given configuration of thinwire piecewise antenna. The numerical procedure 
involved by the algorithm rounds the entries of the Hessian matrix, step by 
step, and the location of the antennas segments can be estimated with high 
precision, so obtaining the image of the shape taken by the antenna. The pa- 
per is structured as follows. In section 2 the mathematical framework about 
the proposed method is presented. In section 3 some numerical results about 
reconstruction examples are discussed, then a conclusion about the method 
results and ongoing work complete the paper. 

2. Numerical approach 

The shape of a biological thin micro-structure can be estimated by means 
of an appropriate embedded emitting antenna and by measuring the spatial 
distribution of the opportune radiated electromagnetic field component. In 
particular, the magnetic field may be the best choice together with the se- 
lection of the appropriate work frequency range. The proposed approach is 
a typical inverse problem: the measured spatial distribution of the magnetic 
field emitted by a VHF thinwire antenna (i.e. with radius less than one hun- 
dredth of the maximum work wavelength) is the input, and the source shape 
reconstruction is the final target. The antenna is assumed to be a piecewise 
linear structure, i.e. a sequence of N^ branches with uniform radius r^^^^ 
and conductivity 



o'omega ■ The antenna is fed at the point Po{xo,yo, Zq) by a sinusoidal current 
source with known frequency/ and amphtude I source- The spatial distribution 
of the magnetic field is detected by a sensors panel with M sensors distributed 
on its surface. As a first approximation, the surrounding medium is supposed 
to be homogeneous, isotropic, with conductivity a, electric permittivity e and 
magnetic permeability fi. Under these assumptions, the shape reconstruction 
involves the minimization of an objective function shown in equation (1), by 
using an iterative procedure. Roughly speaking, step by step the iterative 
procedure looks for a set of segments coordinates that minimize the difference 
between the computed and measured magnetic field at sensor points. 

minp {\\Hcomp{p) - ^means(p) 1 12} (l) 

In equation ([l]) p = [(xi, yi, Zi)...{x]\!^, i/n^, z^^)]'^ is the vector of the carte- 
sian coordinates of the antenna segments ends to be estimated, Hcomp is the 
vector of the magnetic field values at sensors points computed by assuming 
that the segments ends are located at p, and Hmeas is the vector with the 
measured magnetic field values. The values of Hcomp are calculated by the 
forward solver. The forward solver computes the amplitude of the magnetic 
field at a set of points, for a given signal source and antennas characteristics. 
More precisely, the current distribution along the antenna is evaluated by 
solving an appropriate integral equation in frequency domain, derived from 
Maxwells equations. Then, by using the relation that link the currents to the 
magnetic field, this latter can be calculated at the sensors points. The equa- 
tions of the forward model is solved numerically. For the problem addressed 
in this case, the moments method (MoM) via point-matching procedure in 
the frequency domain, is performed . Note that, by this method the integral 
equations involved include the boundary conditions, so that only the emit- 
ter have to be discretized but not the problems domain. This means that 
each antennas branch needs to be split in a finite number of linear segments 
[9]-[Tl]. The problem formulation leads to a first kind Fredholms integral 
equation. In fact, by expressing the electric field E with the retarded mag- 
netic vector potential A, the source frequency {oj = 27r/) and the scalar 
potential 0, as shown in equation (2): 

E = -{\/(l) + juA) (2) 

by introducing the boundary conditions by means of the per-unit-length sur- 
face impedance Zg, the following equation holds: 
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in which tg indicates the component tangential to the wire surface, Is is the 
longitudinal current flowing into the conductor, concentrated on its axis (u' 
as the unit vector) because of the thin- wire assumption [9], [12], u is the unit 
vector tangential to the conductors surface, as shown in figure 1. 
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Figure 1: General thin-wire geometric references. 



By using the relation between the magnetic vector potential A and the 
conduction current density, and the relation between the scalar electric po- 
tential (p aiid the free electric charges density, the well known electric field 
integral equation (EFIE) is obtained. Then, by integrating the EFIE along 
the surface of the conductor, the following modified EFIE is obtained [9], 



j^ II^- In f^^s{l')g{r, r')dl'dl 
+ lLlki-Jn'-^9ir,r')dl'dl 

Zs^^Is{f)dl 



(4) 



Equation (|4]) is a general relation that depends only on thelongitudinal cur- 
rent and on geometrical quantities of the conductors constituting the thin- 
wire structures to be analyzed. In fact, Vt is the length of the exciting con- 
ductor, L is the length of the induced conductor, r and r' are space position 



vectors of the observation and source points, respectively, g{r,r') = ^- r- 

is the Greens function in an unbounded region. The quantity e = e + jua 
take into account the complex medium permittivity and k = \f—uP^e is the 
wave number. As already underlined, the forward problem, represented by 
equation (li]), is numerically solved by splitting the thin-wire antenna into a 
finite number of linear segments of length A. In scientific literature numer- 
ous results |9], [12] show that an acceptable accuracy can be obtained for the 
solution of equation (4), by assuming A < 0.05A and a linear distribution 
of current along each segment. In this way, a linear system of order n + N^ 
is obtained with n as the number of segments. Once the currents are com- 
puted, the magnetic field components in the surrounding medium are given 
by the dipole theory, by superposing the effects of all segments [9], [12]. As 
discussed above, once the direct solver is obtained, the inverse problem can 
be then approached. In figure 2 the blocks diagram of the algorithm used to 
solve the inverse problem is shown. 
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Figure 2: Blocks diagram of the inverse problem algorithm. 

The unknown end-points segments coordinates s are obtained through 
an iterative procedure which, at each step, computes the correction factors 
needed to obtain the new set of coordinates. The correction factors are 
evaluated by the optimization functionality, on the basis of the difference 
in ^-norm between the measured fields and the fields computed by the for- 
ward solver. Since the problem is strictly nonlinear, the solver of the inverse 
problem was based on the Levenberg-Marquardt (LM) algorithm [15]. In 
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brief, the LM algorithm starts from an initial set of reasonable end-points 
segments coordinates s^^^: this set is then updated at each step k by adding 
the solutions As^^^ of the linear system shown in equation A5\ 



J i ~ (^) 

In equation ( 5 ) the matrix B is an approximation of the Hessian matrix 

of F at step k and is equal to G^''^ G^''^ with G^^^ as the gradient of the 
error function F at step k, while /S'-'^^ is an adaptive parameter. Once A s^'^^ 
is known, the new coordinates are calculated by the following equation: 



(fc+i) ^ ^{k) _^ ^g(fc) (g) 



s 



3. Numerical results 

In order to validate the capabilities of the proposed method, some nu- 
merical experiments concerning the position estimation estimation of the 
antennas branches are carried out. The thin-wire antenna is assumed to be 
made by a nickel titanium thin wire 0.1 mm in radius, 2.1 cm in length with 
an electrical conductivity of 1.1 ■ 10^ S/m. In this way the thin wire assump- 
tion is verified. The nickel-titanium alloy is selected because it is the most 
diffused alloy in the biomedical applications. It has to be underlined that 
this type of antenna can be effectively realized in practice and some of the 
human cavities may be investigated. Moreover, the number of branches for 
the antenna is set to N^ = 4, with shape shown in figure 3 and figure 4 by 
the red dashed line with squared markers. The antenna is fed by a sinusoidal 
current source with frequency equal to 100 MHz, and amplitude spacing in 
the following set of values: 0.1 mA, 1 mA, 5 niA, 25 mA and 50 mA. It has to 
be underlined that all these current values as well as the selected frequency 
can be well tolerated by the human tissues. 

The ground surface of the antenna lies on the Y— Z plane as showed in fig- 
ure 3 and figure 4, and the feed point is placed at a distance of 2 mm from 
the ground plane. Two sets of 12 sensors are considered. The sensors are 
assumed to be sensitive to the component along the Z axis of the magnetic 
field and placed on a flat surface in a first case study. 

A semi-cylindrical surface sensors panel is then also considered (figure 4). 
The flat sensor panel is 3.5 cm height and 1.5 cm width, while the semi- 
cylindrical one has 1.5 cm radius and 3.5 cm height. For the above described 
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Figure 3: Conductors position estimated by flat sensor panel: the initial position is that 
a of straight thin- wire antenna 2.1 cm length, perpendicular to x — plane. 
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Figure 4: Conductors position estimated by semi-cylindrical sensor panel: the initial po- 
sition is that of a straight thin- wire antenna 2.1 cm length, perpendicular to x = plane. 



set up, the measured values are calculated by means of the MoM based for- 
ward solver described in section 2. The initial set of end-points segments 
coordinates s^'^^ for the LM optimization algorithm is that of a straight thin- 
wire antenna 2.1 cm length and perpendicular to the ground plane. From 
these coordinates the true position of each branch is then estimated. Fig- 
ure 3 shows the estimated position of the conductor in the case of the fiat 
sensors panel. The picture shows that the distance between the estimated 
position, the black dashed line with plus markers, and the true position, the 
red dashed line with square markers, is close to zero. 

For this configuration. Table 1 shows the ^norm relative error between 
true and estimated coordinates values for a given level of current, and the 
number of iterations. The simulations results show that the error is close 
to zero, also for small values of the source current. This gives benefit for 
the future development of the antennas power system and a benefit for the 



Table 1: Relative error and number of iteration required with the flat sensors panel 
Current [mA] Relative error No. of iterations 



0.1 


6.15e-ll 


169 


1 


6.10e-ll 


165 


5 


7.88e-09 


87 


25 


3.78e-10 


28 


50 


1.53e-ll 


8 



patient. Figure 4 shows the estimated position of the conductor for the semi- 
cyhndrical sensors panel. The picture shows as the distance between the 



Table 2: Relative error and number of iteration required with the semi-cylindrical sensors 
panel 



Current [mA 


Relative error 


No. 


of iterations 


0.1 


3.05e-10 




25 


1 


1.04e-09 




17 


5 


1.75e-ll 




12 


25 


2.55e-12 




8 


50 


l.OOe-09 
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estimated position, the black dashed line with plus markers, and the true 
position, the red dashed line with square markers, is close to zero. For this 
scenario, the ^-norm relative error between true and estimated coordinates 
for a given level of source current, and the number of iterations are shown in 
Table 2. The simulations results show that the error is close to zero also in 
this case. 

Note that a lower relative error has been reached for the flat sensors 
panehthis result can be justified because with the semi-cylindrical sensors 
panel the measured magnetic field components contain more information 
about the real field distribution. Moreover the number of iterations decreases 
as the source current amplitude increases for a same sensors configuration. 
However, the shape of the sensors panel influences also the convergence of the 
optimization algorithm: the semi-cylindrical sensors panel outperforms the 



flat sensors panel, especially for low levels of antenna source currents. This 
result suggests that the number of iterations can be reduced via a proper sen- 
sors panel design, without compromising the precision and without increasing 
the antenna currents, this latter aspect is fundamental for the invasiveness 
of the method. 

4. Conclusions 

In this paper, a numerical method to evaluate the shape of micro-structures 
for bio-medical application with a very low invasiveness for the human body, 
is proposed. A flexible thinwire antenna radiates the VHF waves and then, by 
numerically solving a typical inverse problem, the estimation of the antenna 
location enables to reconstruct the micro-structures image. The typical in- 
verse problem is solved, and first simulation results assess the validity and 
the robustness of the proposed approach. 
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